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Linear analyses for the stability of radial and nonradial 
oscillations of massive stars 
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ABSTRACT 

In order to understand the periodic and semi-periodic variations of luminous O- B- A- 
type stars, linear nonadiabatic stability analyses for radial and nonradial oscillations 
have been performed for massive evolutionary models (8M Q — 90M Q ). In addition to 
radial and nonradial oscillations excited by the kappa-mechanism and strange-mode 
instability, we discuss the importance of low-degree oscillatory convection (nonadia- 
batic g~) modes. Although their kinetic energy is largely conhned to the convection 
zone generated by the Fe opacity peak near 2 x 10 5 K, the amplitude can emerge to 
the photosphere and should be observable in a certain effective temperature range. 
They have periods longer than those of the radial strange modes so that they seem to 
be responsible for some of the long-period microvariations of LBVs (S Dor variables) 
and a Cyg variables. Moreover, monotonously unstable radial modes are found in 
some models whose initial masses are greater than or equal to 60M Q with Z — 0.02. 
The monotonous instability probably corresponds to the presence of an optically thick 
wind. The instability boundary roughly coincides with the Humphreys-Davidson limit. 

Key words: stars: evolution - stars: oscillations - stars: massive - stars: early-type 
- supergiants. 



1 INTRODUCTION 

It is known that various instabilities occur in very lu- 
minous stars. Most luminous stars called Luminous Blue 
Variables (LBVs or S Dor variables) show major events 
(SD(S Dor)-eruptions and SD-phases) with which mass 
loss fro m the star is thought to be enhanced greatly 
(see e.g.. lHumphrevs fe David son 19941 : Ivan" GenderenlliuOll . 
for reviews). Their distribution on the HR diagram 
seems to be bounded by th e Hu mphreys-Davidson (HP ; 
iHumphrevs fe Davidsonll 19791 ) limit. iKiriakidis et all (| 19931 ) 
found that the instability boundary for the radial strange 
modes roughly coincides with the HD limit. They proposed 
that the strange mode instability would yield a strong mass 
loss causing the HD limit. Several nonlin ear analyses for ra- 
dial strange modes have been performed dDorfi fe Gautschvl 
l200d : IChemigovski et aljliool : iGrott et al.ll2005l ) to gener- 
ate pulsation driven mass- loss, but the results seem to re- 
main Tncon£hrave.(In this context, it is interesting to note 
that lAerts et al.l l|2010l ) found a luminous B star to change 
its mass-loss rate on a timescale of the period of photo- 
metric and spectroscopic variation, suggesting a correlation 
between mass-loss rate and pulsation.) 

It is fair to mention that various instabilities other than 
strange mode instability are also proposed for the cause of 
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eruptive mass loss (see e.g., Ivan Gendered lilOll . for a re- 
view). In the present paper ( £I3.2[) we will discuss the pres- 
ence of a monotonously unstable mode in models more lu- 
minous than the HD limit. 

In addition to the major variations (SD-eruptions and 
SD-phases) on timescales of years and decades, LBVs show 
quasi-periodic microvariations on timescales of weeks to 
months. Those variations occur also in non-LBV luminous 
OBA supergiants. They are loosely called a Cyg variables, 
and t he variations are thought to be cau sed by stellar o scilla - 
tions. iDziembowski fe Slawiriskal (|2005l ) and lFadevevI (120101 ) 
claimed these microvariations to be caused by radial strange 
mo de pulsations. 

lLamers et alj l|l99Sl ) argued, however, that periods of 
some of the microvariations in LBVs are much longer than 
the periods of strange modes, and that the positions of some 
of those variables on the HR diagram contradict the strange 
mode instability region. They concluded microvariations in 
LBVs to be consistent with nonradial g-modes. 

In the present paper, we find those longer period vari- 
ations can be interpreted by oscillatory convection (nona- 
diabatic g~) modes associated with the envelope convec- 
tion zone generated by the Fe-opacity peak around T ~ 
2 x 10 5 K. The linear convection modes, which are dynami- 
cally (monotonously) unstable in the adiabatic analysis (g - - 
modes), become overstable (or oscillatory) if nonad i abatic 
effects are included as found bv lShibahashi fe Osakil (|l98ll ) 
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for high-degree (I > 10) modes. These modes have not been 
paid much attention to before, because they are not expected 
to be observed due to the high I values and their amplitude 
being confined to the convective zone. The present paper will 
show some low-degree (£ ^ 2) overstable convection modes 
having periods longer than those of radial strange modes 
can emerge to the stellar surface and hence be observable. 

Before the emerg ence of the OPAL opacity 
l|Rogers fc Iglesiasl |l992), OB-type stars were consid- 
ered to have mostly radiative envelopes. But in reality, the 
Fe-opacity peak at T ~ 2 x 10 5 K in the OPAL opacity is 
strong enough to produce a convective zone of considerable 
thickn ess in the envelopes of OBA stars. ICantiello et al.l 
(2009) argued the importance of the sub-photospheric 
convection in massive stars for various photospheric veloc- 
ity fields including pulsations and large microturbulence 
needed in spect r oscop ic analyses of OB stars. Furthermore, 
iDegroote et al.l l|2010l ) found evidence of solar-like oscilla- 
tions in a O-type star that indicates stochastic excitation 
of oscillations by turbulent convection to work even in hot 
stars. 

There are also less luminous B supergiants which show 
periodic microv ariations with relatively short periods of on e 
to a few days (|Waelkens. et al.lll99Sl ; iLefever et ai]|200it ). 
Some of these variations can be interpreted as supergiant 
extension of SPB stars (SPBsg), which is possible because 
g-mode oscillations being reflected at a convection zone 
above a hydrogen burnin g shell so that strong dissipation 
in the core is supp ressed (|Saio. et al.ll2006l ; lGautschy| [2009; 
iGodart et aLll2009h . 

The present paper discusses the stability of radial and 
nonradial oscillations in massive main-sequence and post- 
main-sequence models, and compares periods of excited 
modes with observed periods of microvariations in super- 
giants. 



2 MODELS AND ASSUMPTIONS 

In order to obtain unperturbed models for stability anal- 
yses covering the upper part of the HR diagram, evolu- 
tionary models in a mass range of 90 ^ Mi/MQ > 8 
were calculated from ZAMS to logT c ff ~ 3.7, using a 
code based on the Henyey m ethod with OPAL opacity ta- 
bles (llglesias fc Roger^ll996l ). where Mi means initial mass. 
The mixing length is assumed to be 1.5 pressure scale 
hei ghts (unles s othe rwise stated). Wind mass loss based 
on IVink et al.l (|200ll ) is included for the models with ini- 
tial mass of Mi 3OA/0. Stellar rotation and core over- 
shooting are disregarded for simplicity. Two sets of chem- 
ical compositions are adopted in systematic calculations; 
(X,Z) = (0.70,0.02) and (0.716,004). An additional model 
sequence was calculated for Mi = 70Mq with the compo- 
sition (X,Z) = (0.71,0.01) to clarify the effect of heavy 
element abundances on the stability of oscillations. 

Table [1] gives stellar parameters at ZAMS, TAMS (the 
end of main-sequence), and at logT c g = 4.0 for selected 
model sequences. It is apparent that significant mass is lost 
during evolution for Mi > 50M Q with Z = 0.02. For the low 
metallicity (Z = 0.004) mass loss is considerably weaker. 

Nonadiabatic radial and nonradial pulsation anal- 
yses were performed using the formulae described in 



Table 1. Parameters for selected models at ZAMS, TAMS and 
at logT Gff =4.0 
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t(10 6 yr) 


2.895 


3.181 


3.747 


4.389 


5.637 


M/Mq 
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logL 
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*(10 6 yr) 


2.925 


3.273 


3.845 


4.556 


6.007 


M/Mq 


79.82 


63.47 


46.68 


37.92 


28.91 


logL 


6.246 


6.100 


5.894 


5.748 


5.548 



ISaio. Winget fc Robinson] (|l983h and lSaio fc Coxl l|l980l ). re- 
spectively. The Lagrangian perturbation of the divergence of 
the convective flux is neglected. The temporal variation of 
radial and nonradial pulsations is expressed as exp(iot) with 
a being eigenfrequency whose imaginary part determines the 
stability;i.e, the pulsation is excited if the imaginary part of 
a is negative. In nu merical ca lculations the frequency is nor- 
malized as lj = <t/*J GM/R 3 with the gravitational constant 
G, stellar mass M, and radius R. 

Although wind mass loss is included in evolutionary 
models for Mi ^ 3OA/0, for oscillation analyses a reflec- 
tive mechanical outer boundary condition (<5P gas — ► 0) is 
employed because theory for the boundary condition in the 
presence of wind has not yet been developed. 



3 EXCITED MODES 

Fig.[T]shows various instability boundaries (visibility bound- 
aries for oscillatory convection modes) obtained in this paper 
and selected evolutionary tracks for a standard composition 
of (X, Z) = (0.7, 0.02). For comparison, positions of periodic 
and semi-periodic supergiant variables are plotted from var- 
ious sources of the observational data. Interestingly, most of 
those variable supergiants reside in the unstable regions. 

The short-dashed line in Fig. [T] shows the instability 
boundary of low-order radial and nonradial p-mode oscilla- 
tions. In the long nearly vertical region of 4.6 > logT e ff > 
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Figure 1. Instability boundaries of various modes with selected 
evolutionary tracks. Black short-dashed and long-dashed lines are 
for low-order and high-oder (radial and nonradial) p-modes, re- 
spectively. Dotted line is the boundary above which monotonously 
unstable modes present. Red and blue solid lines are for nonra- 
dial g-modes with 1=1 and 1 = 2, respectively. Dash-dotted 
lines are for oscillatory convection modes (red line for I = 1, 
blue line for 1 = 2) indicating the boundaries of the region where 
the ratio of the photospheric to the maximum amplitude in the 
interior is larger than 0.2. The number written at the start of 
each evolutionary track indicates the initial mass M; in solar 
mass units. Wind mass loss incl uded for th e evolu tionary mod- 
els with Mi > 30Af Q based on IVink et al.l 1 120011 ). Also shown 
are positions of sup er-giant variable st ars; variable B-type su- 
pergian ts (blue dots [ Lefever et al,l 120071 ). S Dor variables (cyan 
squares |van_j5gnderen| 120011), and a Cyg variables (red triangles 
Ivan Leeuwen. et al.lll99SI : " Ivan Genderen fc Sterke3l2002ft . 



4.3, the kappa-mechanism at the Fe opacity bump ex- 
cites low-order modes, which correspond to fi Cep vari- 
ables llKiriak . dis et alii 19921 : iMoskalik fc Dziembowskll 19921 : 



IPamvatnvkhl 



1900 



The instability boundary bends and becomes horizon- 
tal due to strange mode insta bility l|Kiriakidis et al.lll993l : 
Idatzelll 19941 ; Is aio. et all 199ct ) which occurs when the lumi- 
nosity to mass ratio is sufficiently high (L/M > 10 4 Lq/Mq) 
and radiation pressure is dominant at least locally in the en- 
velope. The insta bility boundary is essen tially the same as 
that obtained by iKiriakidis et all (|l993h for Z = 0.02. In 
the hottest part, the strange modes are associated with the 
Fe opacity peak, while in relatively cooler part the contribu- 
tion from the opacity peak at the second helium ionization 
becomes important. 

The vertical part around log T e s ~ 3.8 is the well known 
blue edge of the Cepheid instability strip. 

A narrow vertical region indicated by a long-dashed line 
around IogT„g ~ 3.8 — 3.9 in Fig. [T] indicates the excitation 
of relatively high order radial and nonradial modes due the 



strange mode instability associated with the hydrogen ion- 
ization zone. We will discuss these modes in i]3.3l 

Monotonously unstable modes exist above the dotted 
line in the most luminous part of Fig. [T] which will be dis- 
cussed in §3.2\ 

In the regions surrounded by blue and red solid lines, 
nonradial g-modes are excited by the Fe opacity bump; red 
and blue lines for 1=1 and I = 2 modes, respectively. Those 
g-modes can be excited even in the post-main-sequence mod- 
els if the mass exceeds 12Mq, because a convection zone 
above the hydrogen burning shell can block some g-modes 
from penetrating to the core. 



3.1 Oscillatory convection modes 

The dash-dotted lines in Fig. [T] are the ranges where oscilla- 
tory convection modes would be observable; red and blue 
lines are for 1=1 and I = 2 modes, respectively. The 
observability is determined based on the ratio between the 
photospheric amplitude and the maximum amplitude in the 
interior, which is discussed below. 

Fig.[2]shows distributions of radial displacements (fr/V, 
bottom panels) and kinetic energy density (det/dr, top pan- 
els) of oscillatory convection modes of £ = 1 in two (hot and 
relatively cool) models from the Mi = 6OM (Z = 0.02) evo- 
lutionary sequence. Two convection modes for each model 
are shown with selected physical quantities. A convection 
zone can be recognized as a zone where the ratio of radia- 
tive to total energy flux is less than unity ;i.e., F ra d/-Ptot < 1- 
The Fe-convection zone occurring in a temperature range of 
5.4 > logT > 4.9 has a considerable thickness but con- 
tains very small mass (< 1O _3 M0 for the hotter model, and 
~ 4x 1O _3 M0 for the cooler model). Because the gas density 
is very low in the envelope of a massive star, the convective 
flux is less than 50% of the total energy flux. In addition, 
the low density makes the gas pressure much smaller than 
the radiation pressure; i.e., (5 C 1 with fi being the ratio of 
gas to total pressure. 

Two modes shown for each model in Fig. [2] are the two 
shortest period convection modes that are most likely ob- 
servable; i.e., ratios of the photospheric amplitude to the 
maximum amplitude in the interior are largest (see also 
Fig. [3| . The numbers in parentheses in the bottom panels 
of Fig. [2] are normalized (by */ GAI/R 3 ) complex eigenfre- 
quencies; the real part represents oscillation frequency, and 
imaginary part growth or damping rate. When the imag- 
inary part is negative the oscillation grows. The eigenfre- 
quencies shown in Fig. [5] indicate that the growth times are 
comparable to the periods. The periods are comparable to 
the periods of g-modes. 

The kinetic energy of the oscillatory convection modes 
is confined to and slightly above the convection zone in the 
hotter model, while it is well confined to the convection zone 
in the cooler model. For the convection modes shown in 
this figure, the oscillation amplitude at the stellar surface is 
comparable to that in the convective zone, indicating those 
modes are very likely observable. 

Since oscillatory convection modes have short growth- 
times, we expect that they have substantial amplitudes in 
the convection zone. For such a mode to be observable, 
however, the amplitude at the stellar surface should be, at 
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Figure 2. Runs of radial displacements (bottom panels) and kinetic energy distributions (top panels) as function of fractional radius 
(r/R) for selected oscillatory convection modes of i = 1 in two supergiant models with different effective temperatures but with the same 
initial mass of 60A/q. Also shown are runs of temperature and the ratios of radiation to total flux (i^rad/^tot) and radiation pressure to 
total pressure (1 — /3). The numbers in parentheses are the real and imaginary parts of eigenfrequencies normalized by 
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Figure 3. Some properties of oscillatory convection modes versus 
effective temperature along the evolutionary track of Mi = 6OA/0 
with Z = 0.02; bottom panel: Real part of normalized frequencies; 
middle panel: growth rate; top panehthe ratio of the amplitude 
of the photosphere to the maximum amplitude. Triangles and 
squares are for 1=1 and I = 2, respectively. Red inverted trian- 
gles are for 1=1 modes in models calculated with a mixing-length 
of one pressure scale-height. 



least, considerable relative to the amplitude in the convec- 
tion zone. We assume that the possibility of detecting these 
convection modes can be measured by the ratio of the pho- 
tospheric amplitude to the maximum amplitude in the in- 
terior. The ratio differs for different modes in a model and 
changes as the stellar parameters change. The top panel of 
Fig. [3] shows variations of the ratio as a function of effective 
temperature along the evolutionary track of Mj = 60Mq 
for a few convection modes of I = 1 (triangles) and I = 2 
(squares) . This figure shows that although many convection 
modes exist in a model, only one or two highest-frequency 
modes (for each £) have large values of the amplitude ra- 
tio and hence are potentially observable. Generally, the am- 
plitude ratio of each mode has a broad peak as a func- 
tion of effective temperature, and the peak value tends to 
be larger for a mode with larger frequency (bottom panel) 
and larger growth rate (middle panel). In other words, the 
visibility is highest for largest frequency modes with high- 
est growth rates in a certain effective temperature range 
(4.2 < logT cfi < 4.6 for the M» = 60M© models). 

Fig. 2] shows contours of the amplitude ratio for the 
most visible convection modes for £ = 1 (solid lines) and 
1 = 2 (dashed lines). Generally, the visibility of the oscilla- 
tory convection modes is better in more luminous and hotter 
models. We assume in this paper that these modes would be 
well visible when the ratio exceeds 0.2. The contours for the 
ratio of 0.2 are shown in Fig. \T\ 

The top panel of Fig. U shows the contours obtained 
for models with (X,Z) = (0.716,0.004). Since the opacity 
is much lower than that for the case of Z = 0.02, the con- 
vection zone around the Fe opacity peak is less extensive. 
Nonetheless, observable oscillatory convection modes with 
low degree I = 1, 2 persist, although the contour for a value 
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Figure 4. Contours of 0.5, 1/3 and 0.2 for the ratio of the am- 
plitude at the photosphere to the maximum amplitude in the 
interior. Solid and dash-dotted lines are for 1=1 and 1 = 2 
modes, respectively. The bottom panel is for normal metallicity 
of Z = 0.02 and the top panel for a metal-poor composition of 
Z = 0.004. 

of the amplitude ratio shifts upward by AlogL ~ 0.8 as 
compared to the case of the standard composition. This in- 
dicates that oscillatory convection modes can cause observ- 
able variations even in the SMC if the stars are luminous 
enough. 

The growth rates of the oscillatory convection modes 
(middle panel, Fig[3]) are very high, especially in relatively 
cooler models; they tend to decrease as T e g increases. The 
reason is probably that in hotter models kinetic energy dis- 
tribution shifts outward and hence a large fraction of the 
kinetic energy resides above the convective zone boundary 
losing the driving. 

It has been thought that the pulsation modes with the 
highest growth rates are as s ociated with stran ge modes (eg. 
iGautschv fc Glatze]|[T990bl ; ISaio. et al.ll 19981 ). The growth 
rates of the convection modes are even larger than those of 
radia l strange modes in ma ssive stars (|Glatzel fc Kiriakidisl 
1 19931 ; IKiriakidis et al.|[r993T ). The oscillation frequencies of 
convection modes lie in the g-mode range as seen in the 
bottom panel of Fig. [3l indicating the periods of convection 
modes are longer than radial modes. 

Fig.[5jshows Q-values of convection modes with the am- 
plitude ratio larger than 0.2 (filled symbols) compared with 
excited radial modes (open circles) for some selected mass 
models, where the Q- value is defined as II ^/p/p with II 
and p being, respective ly, pulsation period and mean den- 
sity (see e.g. lCox|[T980l '). The Q-value of a convection mode 
decreases rapidly as the effective temperature decreases in 
the hottest and the coolest T e a ranges, while it changes lit- 
tle in the intermediate range of effective temperature. At a 
given effective temperature, Q-values of oscillatory convec- 
tion modes tend to be smaller in more massive models. 



Figure 5. Q-values of convection modes (filled triangles for i = 1, 
filled squares for 1 = 2) and excited radial modes (open cir- 
cles) versus effective temperatures along the evolutionary tracks 
of Mi = 25, 40, 60 and 7OM . 

In very massive (M, = 60, 70Mq) models, radial pul- 
sations are excited even in cooler (logT c fi < 4.2) models 
due to mainly the strange mod e effect that works around 
the second He ionization zone (|Glatzel fc Kiriakidisl [l993l ; 
IKiriakidis et al.lll993l ). The Q-values of these strange modes 
increases as the effective temperature decreases; i.e., as the 
depth of the He II ionization zone increases. Fig. [Syndicates 
that periods of oscillatory convection modes are much longer 
than those of radial pulsations in hotter models, although in 
cooler and massive m odels the differences are not very large. 
lLamers et alj |l998) found that the periods of microvaria- 
tions of LBVs are orders of magnitudes longe r than those 
predicted for strange modes by IKiriakidis et al.l (|l993l ). This 
indicates that oscillatory convection modes would be respon- 
sible for long-period microvariations in LBVs (see Sjjbelow). 

Red inverted triangles in Fig. [3] present the property 
of oscillatory convection modes of £ = 1 in models cal- 
culated with a mixing-length of one pressure scale-height 
rather than 1.5 pressure scale-hight adopted in standard 
models. A smaller mixing-length makes the super-adiabatic 
temperature gradient in the Fe-convection zone larger. In 
the models with the smaller mixing- length, real parts of 
eigenfrequencies of oscillatory convection modes tend to be 
smaller, and the amplitude tends to be slightly more con- 
fined to the convection zone. Fortunately, the qualitative 
property of the oscillatory convection modes is insensitive 
to the mixing-length parameter. 

3.2 Monotonously unstable mode 

In very massive models (Mi ^ 60Mq) a monotonously un- 
stable radial mode (with a purely imaginary eigenfrequency) 
is found in the range on the HR diagram indicated by a dot- 
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Figure 6. Growth-time versus effective temperature for 
monotonously unstable modes in models with Mi = 60Mq , 70Mq 
and 9OM . 

ted line in Fig. [T] Fig. [6] shows growth time as a function 
of effective temperature for those unstable modes in models 
with different initial masses. The growth times tend to be 
shorter for more massive stars, ranging from a month to a 
hundred years, which are much faster than the evolution- 
ary change. It is interesting to note that the range of the 
growth times includ es the timescales of S Dor (SD)-phase 
(van Gendered [200lh . To the author's knowledge, the pres- 
ence of monotonously unstable modes in very massive stars 
was not known before. The presence was unrecognized in 
the non-linear radial pulsation analyses performed before 
by various authors, probably because the monotonously un- 
stable mode was concealed by the more fast growing strange 
mode pulsations. 

Amplitude and kinetic energy distribution of the 
monotonously unstable mode in a model with Mi = 70Mq 
are shown in Fig. The amplitude and kinetic energy have 
a peak around the bottom of the convective zone generated 
by the Fe-opacity peak; in this model the convective zone 
ranges from r/R — 0.20 to 0.595. Although, the kinetic en- 
ergy, as well as amplitude, is confined around the bottom of 
the convection zone, the amplitude, after attaining a mini- 
mum at r/R « 0.45, gradually increases toward the surface 
above the convection zone. In a substantial range in this 
model the ratio of the gas to radiation pressure is very small 
(i.e., P < 1). Then, 

^L«(l-/3)«1, (1) 

where L ra( j is radiative luminosity and L^dd the local Ed- 
dington luminosity defined as 

L Edd = 4-KcGMr/K (2) 

with c being the velocity of light and M r the mass within 
the sphere of radius r. This suggests that the monotonously 
unstable mode arises because the radiative luminosity is 
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Figure 7. The amplitude of displacement (j£/»i) and kinetic 
energy distribution (de^/dr) and selected physical quantities as 
functions of the fractional radius (r/R) for a model with Mi = 
70Mq, where re is opacity, and other quantities are the same as 
in Fig. [2] 



very close to the local Eddington luminosity in a substantial 
range of the envelope. 

Fig. [S] shows the distribution of /3 in selected mod- 
els, where solid lines are used for models which have a 
monotonously unstable mode, while dashed lines are for 
models without it. There are two minima of f3 at logT « 5.4 
(Fe opacity peak) and at logT ~ 4.6 (He II ionization). Ap- 
parently, /3 must be well below 0.1 in a range including the 
both opacity peaks for a monotonously unstable mode to 
appear. Note that the hotter model of 70Mq (black dashed 
line) has a distribution of /3 around 5.2 < logT < 5.4 similar 
to that of the cooler model (black solid line), but f3 in the 
range 4.5 < logT < 5 is significantly higher than that in 
the cooler model, which inhibits a monotonously unstable 
mode in the hotter model. With a low-metal abundance of 
Z — 0.004 no monotonously unstable mode is found even in 
models of Mi = 9OM0. The reason is obvious from the blue 
dashed line in Fig. [5] /3 at the Fe opacity peak is substan- 
tially larger, although a further increase in mass seems to 
cause monotonously unstable modes. 

The presence of such a monotonously unstable mode 
probably corresponds to the presence of a n op tically 
thick wind as discu ssed by iKato fe Ibenl |l992) and 
iNugis fe Lamersl (|2002h for WR stars. It is known that opti- 
cally t hick winds also occur in nova models. IKato fe Hachisul 
(2009) have shown that for non-extreme nova cases both 
static models and models with winds are possible. The pres- 
ence of a monotonously unstable mode in a massive evolu- 
tionary model might indicate that the static model transits 
to a model with an optically thick wind. 

It is interesting to note that the boundary in the 
HR diagram for the presence of monotonously unstable 
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Figure 8. Distribution of the ratio P ga s/f (= 0) in the interior 
of selected models. Solid lines are used for the models which has 
a monotonously unstable mode, while broken lines for the models 
having none. 



modes (dotted line in Fig. QJ ro ughly coincides with the 
Hum phreys- Davidson (HD) limit l|Humphrevs fe Davidson! 
Il979h : this may suggest the HD limit to be related to the 
presence of optically thick winds. 



3.3 High-oder modes excited at H ionization zone 

In the narrow vertical region around logT c g ~ 3.9 seen in 
Fig. [JJ relatively high order (n ~ 6 — 10) radial as well as 
nonradial modes are excited at the hydrogen ionization zone 
at log T ~ 4. Since the normalized frequency of the excited 
mode increases as the effective temperature decreases, the 
red-edge is formed where the frequency exceeds the criti- 
cal frequency given as ^ \JT\R/ H p with the pressure scale 
height Hp and Ti = (din P/d In p) ac i. (If the frequency ex- 
ceeds the critical frequency, the pulsation is not reflected at 
the outer boundary and is expected to be dissipated in the 
outermost layers.) 

Fig. [9] shows the amplitude distribution and the work 
curve for a high-order radial mode excited in a model rela- 
tively close to the blue edge. Obviously the excitation occurs 
at logT « 4, in the hydrogen ionization zone. The ampli- 
tude of the mode is extremely confined to the surface layers, 
where the radiation pressure is much larger than the gas 
pressure. Since unstable modes of this t ype are present even 
in the non-adiabatic reversible (NAR; Gautschv fc Glatzell 
Il990bl) approximation where thermal-time is artificially set 
to be zero, we can identify such a mode as a strange mode 
trapped above the opacity peak of the hydrogen ionization. 
A similar mode with a similar frequency is excited for 1=1 
and 1 = 2. These excited nonradial mod es probably corre - 
spond to the trapped modes discus sed in IGautschvl (|2009h . 
In addition to the trapped modes IGautschvl ( 20091 ) found 



Figure 9. Amplitude distribution and work curve (top panel) 
for an excited 9th-order radial pulsation. Some physical quantities 
are shown in the bottom panel. 



several untrapped modes to be excited around the trapped 
mode freq uency, by using the Riccati- shooting method dis- 
cussed in iGautschv fc Glatzel (1990a). The present study 
using a finite difference method, h owever, could no t find 
those untrapped modes. According to lGautschvl (|2009l ) those 
modes have many nodes with considerable amplitude in the 
deep interior around r/R ~ 10~ 2 . The finite difference code 
even with ~ 10 4 grid points seems ineffective to find such 
untrapped modes. 



4 COMPARISON WITH OBSERVED PERIODS 

Fig. [10] compares periods of radial and nonradial modes ex- 
cited in 16Mq and 30Mq models with observed periods of 
relatively less luminous {log L/Lq < 5.5) supergiant vari- 
ables as function of T e a. Those v ariables consis t of m ost of 
the periodic B stars analyzed bv lLefever et al.l ((2007) , and 
relatively cooler a-Cygni variables including a Cyg itself. 
Most of the periods of the hotter (logT c fi > 4.15) less lumi- 
nous supergiants seem consistent with oscillatory convection 
modes, SPB-type g- modes (SPBsg), or p-modes. 

Relatively cool (logT e ff ;$ 3.95) stars in Fig. \W\ are less 
luminous a-Cygni variables. They lie in the T e g range where 
high-order p-modes are excited (Fig. [T}, although most of 
the periods are longer than those of high-order p-modes. 
Probably they are pulsati ng in g-modes ex cited in a similar 
range of T e g as found by IGautschvl (|2009h . who found that 
g-modes having periods of ranging from about 10 to 20 days 
are excited at the hydrogen ionization zone, although a large 
fraction of their kinetic energy resides in the deep interior. 

No modes can be assigned to the three stars in a range 
of 4.15 > logT off > 4.05 in Fig. HTJl which ar e the coolest 
members in the data from iLefever et ail l|2007l ) (FigUJ. The 
cause of the oscillations in these stars is not clear. 
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Figure 10. Periods of radial and nonradial modes excited 
in 16 and 3OA/0 models are compared with observed peri- 
ods of relatively less luminous (log L/Lq ^ 5.5 ) supergiants 
(large dots) which consists of periodic B stars from lLefever et al.l 
(2007) and relatively cool er les s luminous a Cyg variables from 
Ivan Leeuwen. et al. I GUI) and Ivan Genderen fc Sterkenl J2002I) . 
The periods of a Cyg itself are adopted from | Luc^"Il9761) ■ while 
the effective temperature from Sch iller fc Przvbillal J2008h . Tri- 
angles and squares correspond to i = 1 and 1 = 1, respectively. 
Filled triangle and squares are convection modes for which the 
amplitude ratios (see A3. Ill are larger than 0.2. Open circles are 
radial modes. Swarms of t = 1,2 modes in a range of 4.35 > 
logT efI > 4.15 for 16M and in a range of 4.5 > logT cff > 4.35 
for 30AfQ are high-order g-modes (i.e., SPBsg-type oscillations). 



Fig. 1 1 1 1 compares periods of radial and nonradial modes 
excited in higher mass models of Mi = 60 and 70M© with 
observed ones for microvariations in luminous (log L/L ^ 
5.5) supergiants which consist of LBVs and luminous a 
Cyg variables. In these massive stars, radial and nonradial 
strange modes are excited (open symbols), and oscillatory 
convection modes (filled triangles and squares) are likely ob- 
servable (the amplitude ratio is larger than 0.2) in a wide 
range of T e g, while SPB-type g-modes are not excited. 

Fig. [TT] indicates that most of the periods of microvari- 
ations of LBVs and luminous a Cyg variables are consistent 
with the periods of oscillatory convection modes and strange 
modes. 

The comparisons between observed and theoretical pe- 
riods in Figs. [10] and [TT] indicate that a-Cygni variables are 
inhomogeneous. Luminous ones are similar to the microvari- 
ations of LBVs which are identified as oscillatory convection 
modes or strange modes, while less luminous cooler ones, 
including a Cyg itself, seem to be g-modes (and possibly in- 
cluding high-order strange modes) excited at the hydrogen 
ionization zone. 



Figure 11. Periods of radial and nonradial modes excited in 
high mass models of Mi = 60M Q and 70M Q with Z = 0.02 are 
compared with observed periods of luminous (log L/ Lq ^ 5.5) 
super giants (red do t s) wh ich consists of microvariations of LBVs 
from [ Lamers et al. I <199Sft and luminous a Cyg variabl e s from 
Ivan Leeuwen. et al.l l ll99a) and Ivan Genderen fc Sterkenl <2002ft . 
Symbols for theoretical periods have the same meanings as in 
Fig. [10] 



M 

° 1 




Figure 12. Periods of radial and nonradial oscillations excited 
in models with different heavy element abundances are compared. 
The symbols for the excited modes are the same as in Fig. 1101 
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Figure 13. The same as Fig. [T]but for a metal-poor composition 
of (X,Z) = (0 . 716,0 04). Also plotted are s ome S Dor variables 
(van Gendered [20011) and a Cyg variables llvan Leeuwen. et al.l 
1 19981) in SMC 

5 LOW METALLICITY CASES 

Since there are some LBVs and a Cyg variables in SMC, it 
would be useful to discuss the results for low metal abun- 
dances. Fig. [T5] compares periods of radial and nonradial 
modes excited in models with Mi = 70M© having differ- 
ent metal abundances. The difference between the cases of 
Z = 0.02 and Z = 0.01 is small, while the case of Z = 0.004 
is significantly different from the former cases. Because of 
the Fe-opacity peak being reduced, strange mode instabil- 
ity for low-order radial and nonradial modes are reduced 
significantly for the Z = 0.004 case, while the excitation 
of SPBsg-type g-modes seems to be enhanced in the lowest 
metallicity case (the reason is not clear). 

It is interesting to note that even for the low metallicity 
of Z = 0.004 oscillatory convection modes still seem to be 
observable in a considerable range of T e g (see also Fig. 2} . 

Fig. [13] is the same as Fig. [l] but for Z = 0.004; i.e., 
shows instability ranges for radial and nonradial modes and 
observable ranges for oscillatory convection modes in the HR 
diagram. Compared to the case of Z — 0.02 (Fig.[T]), the in- 
stability region for low-order radial and nonradial p-modes 
(i.e., P Cep instability strip) has disappeared completely. 
The instability boundary for the low-order strange mode in- 
stability has shifted upward conside rably, which agrees with 
the result of Kiriakidis ct al. | (|l993h . 

On the other hand, less affected are the SPBsg-type g- 
mode instability regions and the observable region of oscilla- 
tory convection modes. The instability region for high-order 
modes in 3.95 > logT e fF > 3.8 is not affected because they 
are excited at hydrogen ionization zone. Positions of some 
LBVs (S Dor variables) and a Cyg variables in SMC are 
also plotted in Fig. 1131 Fig. [14] compares quasi-periods of 
the microvariations of those stars with theoretical periods 



Figure 14. Periods of excited modes in metal-poor models 
(Z = 0.004) of Mi = 30, 50, 70, and 90M Q as function of T cB 
are compared with observed periods of microvariations of LBVs 
and a Cyg variables in SMC. Symbols for theoretical periods have 
the same meanings as in Fig. 1101 

predicted for metal poor models with initial masses ranging 
from Mi = 30 to 9OM . 

In contrast to the Galactic and LMC cases, oscillatory 
convection modes seem to be responsible for none of the 
SMC stars shown in these figures. Fig. 1141 seems to indicate 
that periods of low-order strange modes are comparable with 
those of the microvariations of R40 and R45 (and probably 
R42). For the strange mode interpretation to be true, the 
initial mass must be greater than ~ 8OA/0. Fig. 1131 however, 
indicates the initial mass of these stars to be 40 — 50M©; 
i.e., the luminosities of these stars seem too low at least by 
a factor of two. Therefore, the mode identification for the 
microvariations of R40, R42 and R45 is unclear. 

For S18, an a Cyg variable in the SMC, only 
the effective temperature is indicated in Fig. [2] be- 
cause the quasi-periods of its variability are uncertain, 
van Genderen fc Sterkenl (|2002h found three types of vari- 
ations on timescales of a few years, ~ 150 days, and a 
few days. Figs. \T3\ and \T%\ indicate the shortest timescale 
variation to correspond to the SPBsg type g-modes. For 
the longer-periods variations, however, no excited oscillation 
modes can be assigned. 



6 CONCLUSIONS 

The stability of radial and nonradial modes in massive stars 
(90 ^ Mi /Mo ^ 8) was investigated. The K-mechanism as- 
sociated with the Fe-opacity bump at T ~ 2 x 10 5 K excites 
low-order radial and nonradial p-modes in main-sequence 
and post- main- sequence models (/? Cep instability strip) and 
low-degree high radial order g-modes in post-main sequence 
models (supergiant SPB stars; SPBsg ). Interestingly, it is 
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found that under a low-metal composition of Z = 0.004, 
the SPBsg instability region still remains, while the j3 Cep 
instability strip disappears. 

In luminous models with L/M > 1O 4 L©/M (for Z = 
0.02) strange mode instability occurs associated with opac- 
ity peaks due to Fe and He II ionizations. For the metal-poor 
composition Z — 0.004 strange mode instability appears 
only in more luminous models. 

In a narrow strip of 3.95 > logT e fi > 3.85, high-order 
radial and nonradial pulsations are excited at the hydrogen- 
ionization zone. The pulsation amplitude is strongly con- 
fined to the outermost layers. 

Low-£ oscillatory convection modes associated with the 
Fe convection zones in massive stars are found. They are 
likely observable in sufficiently luminous models in a con- 
siderable range of effective temperature. The oscillatory con- 
vection modes tend to have growth times comparable with 
the periods. Their periods are longer than those of strange 
modes or SPBsg-type g-modes, and are comparable to those 
of long-period microvariations of LBVs and luminous a Cyg 
variables. Rapid growth rates of the oscillatory convection 
modes are consistent with the presence of irregularities in 
the microvariations. 

In addition to the oscillation modes, a monotonously 
unstable mode is found in luminous models with logL > 5.9 
and logTefi < 4.3. The boundary in the HR diagram for 
the presence of the monotonously unstable mode roughly 
coincides with the HD limit. The growth rate is more than 
hundred times faster than the stellar evolution rate. It prob- 
ably suggests the presence of an optically thick wind in such 
a luminous star. 

Finally, it should be noted that in the present stability 
analysis the effect of stellar winds is not included in the 
outer boundary condition. In the future, we have to clarify 
how the stellar wind affects the stability of oscillations of 
massive stars. 



ACKNOWLEDGMENTS 

I am grateful to Alfred Gautschy for making sample calcula- 
tions with an alternative numerical approach to confirm the 
presence of overstable low-degree convection modes, and for 
helpful comments on a draft of this paper. I am also grateful 
to the anonymous referee for useful comments, which have 
improved the paper considerably. 



Aerts C, Noels A., Rainer M., et al., 2010, A&A, 519, 
A38 

Dorfi E.A., Gautschy A., 2000, ApJ, 545, 982 

Dziembowski W.A., Slawihska J., 2005, AcA, 55, 195 

Fadeyev Yu.A., 2010, Ast.Let., 36, 362 

Gautschy A., 2009, A&A, 498, 273 

Gautschy A., Glatzel 1990a, MNRAS, 245, 154 

Gautschy A., Glatzel 1990b, MNRAS, 245, 597 

Glatzel W., 1994, MNRAS, 271, 66 

Glatzel W., Kiriakidis M., 1993, MNRAS, 263, 375 

Godart M., Noels A., Dupret M.-A., Lebreton Y. 2009, 

MNRAS, 396, 1833 
Grott M., Chernigovski S., Glatzel W., 2005, MNRAS, 360, 

1532 

Humphreys R.M., Davidson K., 1979, ApJ, 232, 409 
Humphreys R.M., Davidson K., 1994, PASP, 106, 704 
Iglesias C.A., Rogers R.J., 1996, ApJ, 464, 943 
Kato M., Iben I.Jr, 1992, ApJ, 394, 305 
Kato M., Hachisu I., 2009, ApJ, 699, 1293 
Kiriakidis M., El Eid M.F., Glatzel W., 1992, MNRAS, 255, 
lp 

Kiriakidis M., Fricke K.J., Glatzel W., 1993, MNRAS, 264, 
50 

Lamers H.J.G.L.M., Bastiaanse M.V., Aerts C, Spoon 

H.W.W., 1998, A&A, 335, 605 
Lefever K., Puis J., Aerts C, 2007, A&A, 463, 1093 
Lucy L.B., 1976, ApJ, 206, 499 

Moskalik P., Dziembowski W.A., 1992 A&A, 256, L5 

Nugis T., Lamers H.J.G.L.M., 2002, A&A, 389, 162 

Pamyatnykh A.A., 1999, Acta Astron, 49, 119 

Rogers, F.J., Iglesias, C.A., 1992, ApJS, 79, 507 

Saio H., Cox J.P., 1980, ApJ, 236, 549 

Saio H., Baker N.H., Gautschy A., 1998, MNRAS, 294, 622 

Saio H., Kuschnig R., Gautschy A., Cameron C, Walker 

G.A.H., Matthews J.M., et al., 2006, ApJ, 650, 1111 
Saio H., Winget D.E., Robinson E.L., 1983, ApJ, 265, 982 
Shibahashi H., Osaki Y., 1981, PASJ, 33, 427 
Schiller F., Przybilla N., 2008, A&A, 479, 849 
van Genderen A.M., 2001, A&A, 366, 508 
van Genderen A.M., Sterken, C, 2002, A&A, 386, 926 van 

Genderen A.M., Sterken C, Jones A.F., 2004, A&A, 419, 

667 

van Leeuwen F., van Genderen A.M., Zegelaar, I., 1998, 

A&AS, 128, 117 
Vink J.S., de Koter A., Lamers H.J.G.L.M., 2001, A&A, 

369, 574 

Waelkens C, Aerts C, Kestens E., Grenon M., Eyer L., 
1998, A&A, 330, 215 



REFERENCES 

Aerts C, Lefever K., Baglin A., Degroote P., Oreiro R., 

Vuckovic M., Smolders K., Acke B., Verhoelst T., Desmet 

M., et al., 2010, A&A, 513, Lll 
Cantiello M., Langer N., Brott I., de Koter A., Shore S.N., 

Vink J.S., Voegler A., Lennon D.J., Yoon S.-C, 2009, 

A&A, 499, 279 
Chernigovski S., Grott M., Glatzel W., 2004, MNRAS, 348, 

192 

Cox J. P., 1980 Theory of Stellar Pulsation (Princeton Univ. 
Press) 

Degroote P., Briquet M., Auvergne M., Simon-Diaz S., 



